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Optimal Actuator and Observation Location for Time-Varying 
Systems on a Finite-Time Horizon * 
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Abstract 

The choice of the location of controllers and observations is of great importance 
for designing control systems and improving the estimations in various practical prob¬ 
lems. For time-varying systems in Hilbert spaces, the existence and convergence of the 
optimal location based on linear-quadratic control on a finite-time horizon is studied. 

The optimal location of observations for improving the estimation of the state at the 
final time, based on Kalman filter, is considered as the dual problem to the LQ optimal 
problem of the control locations. Further, the existence and convergence of optimal 
locations of observations for improving the estimation at the initial time, based on 
Kalman smoother is discussed. The obtained results are applied to a linear advection- 
diffusion model. 

Keywords: Approximation, Kalman filter, Kalman smoother, linear-quadratic control, op¬ 
timal observation location 

1 Introduction 

The choice of the locations of control hardware, such as sensors and actuators, plays an im¬ 
portant role in the designs of control systems for many physical and engineering problems. 
Proper locations of sensors and actuators is essential to improve the performance of the con¬ 
trolled system. Many researchers have focused on the study of finding the optimal locations 
of control hardware and different criteria of optimising control locations were established, 
such as maximization of observability and controllability fl7l . EH, minimizing the linear 
quadratic regulator cost Ii23l . Geromel Ifl4l successfully reformulated the LQ cost function 
into a convex optimization problem by mapping the locations of controller into zero-one 
vectors and expressed the solution of classic LQ problem in terms of a Riccati equation. 
Morris ll22il optimized controller locations of time-invariant systems on an infinite-time 
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horizon in Hilbert spaces by solving an algebraic Riccati equation and showed the con¬ 
vergence of optimal controller locations of a sequence of approximated finite-dimensional 
systems. Further, an algorithm lITOll for the linear quadratic optimal problem of controller 
locations based on the convexity shown in Ifl4l are introduced. 

The issue of observations is also of great importance of many estimation problems for 
stochastic systems, such as weather forcasting and data assimilation in meteorology. For 
this kind of problems, observations always have low temporal and spatial density. The lack 
of observations is a major bander of preventing the improvement of estimations and leading 
to the accuracy of predictions. On one hand, the insufficient observations become the main 
reason that many works are introduced to improve approaches of estimations in in the recent 
years. On the other hand, one possibility to improve the predictive or estimation skill for 
specific problems is to target the locations of observations which can potentially result in 
the largest forecast improvement in order to make observations more efficient. The better 
choice of locations of observations can help making more progress of the predictive or esti¬ 
mation skills. In contrast, improper observations probably make no sense to the accuracy of 
predictions and lead to the waste of resources by optimizing the improper parameters. There 
are several papers focusing on this problem from the perspective of applications. For finite¬ 
dimensional systems in practice, approaches based on singular value decomposition (Oil. 
HI ) always help determining the direction with the strong influence of observations. How¬ 
ever, it cannot solve the optimal problem of observation locations. Motivated by problems 
of data assimilation in meteorology, we estimate unknown random variables by Kalman 
filter and smoother, which has been theoretical foundation of one of the most popular data 
assimilation approaches in last decade. In fact, since 1960’s, besides of applications in me¬ 
teorology, the Kalman filter and smoother |fl9ll were widely applied in many other fields 
to produce optimal linear estimations of states and parameters through a series of observa¬ 
tions over time. It provides us an opportunity to define and search for optimal locations of 
observations by minimizing the covariance based on Kalman filter and smoother. 

In this paper, we will start from the infinite-dimensional state space to consider the 
optimal location problem of controllers and observations for time-varying systems on a 
finite-time horizon. First, we study the linear-quadratic optimal location control problem 
for both deterministic and stochastic systems and develop conditions guaranteeing the ex¬ 
istence of optimal locations of linear quadratic control problems in Section [2] Associated 
with practical applications, since optimal control problems cannot be solved directly in 
infinite-dimensional spaces, a sequence of approximations of the original time-varying sys¬ 
tem have to be considered. Thus, in Section [3j analogical to the approximation theory 
of time-invariant systems, we introduce the similar approximation conditions of evolution 
operators so as to ensure that the approximated control problems converge to the optimal 
control problem of the original infinite-dimensional time-varying system. Further, we show 
the convergence of minimal costs and optimal locations of the sequence of approximations. 
In Section 0] and Section [5] we derive the Kalman filter and smoother of time-varying sys¬ 
tems in the integral form on Hilbert spaces. Then, by duality between Kalman filter and 
linear-quadratic optimal control, under certain conditions, the nuclearity of the covariance 
can be guaranteed. In Section[6]based on Kalman filter and smoother, the existence and con¬ 
vergence of optimal location observations of the estimation of the model state for stochastic 
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systems is shown. Finally, we apply the obtained results to a three-dimensional advection- 
diffusion model with the special construction of the emission rate in Section [7] In this 
example, the operator splitting technique with spatial and temporal discretization is applied 
to simulate the practical application in meteorology. 

2 Existence of optimal actuator locations 

Throughout this paper, we will always assume that the state space of the time-varying sys¬ 
tem is a real separable Hilbert space X, and the input and output space are Hilbert spaces 
denoted by U and Y, respectively. First, we introduce the notion of mild evolution operators 
for the time-varying system. 

Definition 2.1. Denote : {(f, s)\ — oo<a^s^t^b< oo}. We call T(-, •) : —)■ 

C{X) a mild evolution operator if 

1. T(t, t) = I, 

2. T(t , r)T(r, s ) = T(t, s), a^s^r^t^b, 

3. T(-, s) : [s, 6] —>■ C(X) and T(t, •) : [a, t] —> C{X) are strongly continuous. 

4. A := sup ||T(f, s)|| < oo. 

(M)er£ 

In the following we assume that T(-, •) : T h n —*■ CiX) is a mild evolution operator, and 

B £ if (a, b ; U, X) with B* € Lf (a, 6; X, U). Here 

L£°(a, 6; X,Y) := {F : [a, b] —y C(X. Y) \ F is strongly measurable and 

||F||oo := ess sup ||.F(f)|| < oo}. 

t£[a,b\ 

For an initial time to F [a, 6], we consider the time-varying system described by 

x(t) = T(t,t 0 )xo+ [ T(t,s)B(s)u(s)ds , t€[t 0 ,b], (1) 

Jt 0 

where xq £ X and u £ L 2 (to,b;U). We are intersted in the following linear-quadratic 
optimal control problem. 

Linear-Quadratic Optimal Control Problem: Find for xo £ X a control uq £ 1 2 (to - 6; U) 
which minimizes the cost functional 

f b 

J(t 0 ,x 0 ,u) = (x(b),Gx(b)) + / \\C(s)x(s)\\ 2 + (u(s),F(s)u(s))ds, (2) 

Jt 0 

where the function x is given by dTJ. Here C £ L ( f > {a,b]X,Y), G £ C{X) and F £ 
L f (a, b: U, U ), Fit) is self-adjoint and nonnegative for fixed f, and F- 1 £ L?(a,b-,U,U). 
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It is well known, see lfl5ll . that the linear-quadratic optimal control problem possesses 
for xo £ X a unique solution uq, which is given by uo(t) = —L(t)x(t), t £ [to,b], 
L(t) = F _1 (t).B*(f)n(t), such that the minimum of the cost functional is given by 

min J(t 0 ,x 0 ,u) = J(t 0 ,x 0 ,u 0 ) = (x 0 , n(f 0 )x 0 ), 

ueL 2 (t. 0 ,b;U) 

where the self-adjoint nonnegative operator 11(f) is the unique solution of the first integral 
Riccati equation (IRE) 

n(f)s = T*(b, t)GT(b, t)x 

rb 

+ j T*(s,t)[C*(s)C(s)-U(s)B(s)F-\s)B*(s)U(s)]T(s,t)xds (3) 

and the second IRE 

U(t)x =Tfr(b,t)GT u (b,t)x (4) 

rb 

+ J T^{s,t)[C*(s)C{s) + U(s)B(s)F- 1 (s)B*{s)U{s)}T n (s,t)xds, 

where 

Tn(f, t)x = T(t, t)x — J T(t, s)B(s)F~ 1 (s)B*(s)H(s)Tn(s , r)xds, (f, r) £ T*. 

Now we consider the situation having the opportunity to choose m locations to control 
and each location varies over a compact set Q C K(. We indicate these m locations by 
the parameter r £ Q m , and denote the location-dependent input operator B(-) by B r (■). 
Throughout the rest of the paper, by a time-varying system with location-dependent input 
operator the time-varying system <[T|) and the cost functional © with B r instead of B is 
meant. The corresponding solution of the IRE and the Riccati operator L are denoted by IT,, 
and L r , respectively. 

In most cases, the initial state xo is not fixed. This indicates several different ways to 
define the optimal actuator location problem. We take two possible ways into account here. 
The first one is to minimize the cost with the worst choice of initial value, which is 

max min J r (xn,u) = max (xn, n r .(fn).To) = ||II r (fo)||. 
\\x 0 \\=lu£L 2 (t 0 ,b;U) ||®o||=l 

Let T r (to) := ||n r (t 0 )||, the optimal performance of r is T(to) = inf rS Qm l|n r (f 0 )l|. 

The second one is to assume that the system is stochastic. Thus, we need to consider the 
trace of n r (fo) instead, since the trace indicates the sum of the diviation of the state vector 
in each coordinate. Thus the evaluation of the particular performance of r is given by the 
nuclear norm of n r (fo), which is (fo) = ||II r (fo)||i. Further, the optimal performance is 

h{to) = inf ||n r (f 0 )||i. 

For time-invariant problems on an infinite time horizon this problem was studied in If22ll . 
In this section we prove the existence of optimal control locations for deterministic as well 
as stochastic time varying systems on a finite-time horizon. 
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Theorem 2.2. Let {B r } r , r G Tl m , be a family of compact operator valued functions with 
the property that lim r _» ro \\B r — -B ro ||oo = 0 for some ro G Tl m . Then the solutions of the 
corresponding integral Riccati equations II r satisfy 


lim ||II r (t) — n ro (i)|| = 0, tG [a, b\, 

r —>ro 

and there exists an optimal location r such that for any initial time to G [a, b], 

i(to) = ||n,(to)|| = inf ||n r (to)||. 

Proof Thanks to the assumptions on B r , there exists 5 > 0 such that A b ■= sup{||-B r (f)|| j 
t G [a,b\, ||r — ro|| ^ 5} < oo. We denote by B(ro,<5) the set B(ro,5) := {r G f l m : 
|r — ro|| ^ <5}. Thus, fl5l Theorem 5.1] implies for every x G X 

U r (t)x -»• n ro (t)x, r —>■ ro- 

For any feedback control u(t) = L(t)x(t), L G L£°(a, b; X, U ), 


(x(t), U r (t)x(t)) ^ J(t,x(t),u ) 

r b _ 

= (x{b), Gx(b)) + / \\C(s)x(s)\\ 2 + (L(s)x(s),F(s)L(s)x(s))ds 

= \\G^Tp r (b, t)x(t )|| 2 

+ J \\C(s)T L r (s,t,)x{t )\\ 2 + \\F^ (s)L{s)T Lr (s,t)x(t)\\ 2 ds, 

where T Lr (t, t)x = T(t, r)x + f* T(t, s)B r (s)L(s)T Lr (s, r)xds, ( t , r) G T b a . 

Since the family B r is uniformly bounded by A b on B(ro, 5), ffl5l Theorem 2.1] implies 
for all r G B(ro, S), (t, r) G T^, ||Ti(t,r)|| < A exp(AAfi ||T|| 00 (t — r)). Further, because 
C G L“(o, 6; X, Y), F G Lf 3 (a, 6; U, U ), there exists a constant An, independent of t and 
r G B(ro, 5), such that HF^-Hoo ^ An- 

For S r = C*C — L*FL r , where L r = F~ l B*Ti r , we obtain 

U r (t)x -U ro (t)x = J T*(s,t) (S r (s) — S ro (s))T(s,t)xds, x G X. 

Since F~ 1 G Ff(a. 6; U, U) and the operator B ro (t) is compact for any t G [a, b], we have 

\\L* r (t) - L* ro (t)\\ ^ ||i ?_1 || 0O (||n r (t)||||B r .(t) -B ro (t)\\ 

+ ||(n r (t) - n ro (t))B ro (t)||) —> 0, r — >■ r 0 , 

which shows 

\\S r (t) - S ro (t)\\ < \\L* 0 (t) — L*(t)\\\\F(t)L ro (t)\\ 

+ ||L*(t)*’(t)||||L ro (t) - L r (t )|| —> 0, r —► r 0 . 
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From the uniform boundedness of F, B r and II r on B(ro,<5), L r and further S r are uni¬ 
formly bounded for all t G [a, b] and B(ro, 5). According, thanks to the dominated conver¬ 
gence theorem, we obtain ||n r (t) — II ro (t)|| —>• 0, r —>• Tp. 

Additionally, since r G Q m , <T m is a compact set, there exists an optimal location f such 
that ||n^(t 0 )|| = inf re nm ||n r (t 0 )||. □ 

Theorem 12.21 shows the continuity of optimal actuator locations and existence of the 
optimal location in the operator norm. For stochastic systems, the above problem leads 
to the nuclear norm. Thus, first we develop conditions which guarantee that the Riccati 
operator is a nuclear operator. Similar to (8] Theorem 3.1], we have 

Theorem 2.3. Let T(-, •) be a mild evolution operator on X, B € L“(a, 6; C p , X), and 
C G L“(a, 6; X, C q ). Then for any to G [a, b] we have: 

1. The observability operator Ct 0 : X —> if (to, 6; C q ) defined by 

(Ct 0 xo)(-) = C(-)T(-,t 0 )x 0 , x 0 G X, 
is a Hilbert-Schmidt operator; 

2. The controllability operator Bt 0 : L 2 (fio, 6; O') —> X defined by 

Bt 0 u = I T(b, s)B(s)u(s)ds 
Jt 0 

is a Hilbert-Schmidt operator; 

3. C* Q Ct 0 and B t(t Bf are nuclear operators. 

Proof. 1. Defining C to y : X —> L 2 (fo, b), i G {1,..., q} 


(Ct 0 ,i x o)( s ) = ( C(s)T(s,t 0 )xo,ei), s ^ t 0 , 

where {e;} is the standard orthogonal basis of C q . We have 

|(Cf 0 ,Oo)(»| = |(C , (s)T(s,to)xo,e i )| ^ ||<7(s)T(s,f 0 )z 0 |||N| 

< ll c '( s )llll r (Mo)||||0)|| < oo. 

E51 Theorem 6.12] implies that Ct 0 y is Hilbert-Schmidt, that is, for any orthogonal basis 

{ej of X, we have £? =1 £“ i IlCto,^ ||| 2(t0)6) < oo. Since 

q 

WCtotjWhtofi) = XI H^oG®jlll 2 (t 0 ,fe)’ 

2=1 

we have 

oo q oo 

XI ll^0^illL 2 (t 0 ,b;C9) = XI XI H^ll! 2 (* 0 , b ) < 00 ’ 

j=l i=lj=l 
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which shows that Ct 0 is a Hilbert-Schmidt operator. 

2. According to f26l Theorem 6.9], Bt 0 is Hilbert-Schmidt if and only if B* 0 is Hilbert- 
Schmidt. An easy calculation shows B* 0 : X —? L 2 (to , 6; U), 

From part 1, B* 0 is Hilbert-Schmidt, and so is Bt 0 . 

3. Since ||C t *C to ||i < ||C* 0 \\Hs\\C to ||j?s < oo and \\B^B to \\i < \\B* 0 \\Hs\\J3t 0 \\HS < 

oo, C* C tQ and Bt 0 B* 0 are nuclear operator. □ 

Corollary 2.4. Assume that the input space U and the output space Y are finite-dimensional 
and G that is a nuclear operator, then the unique nonnegative self-adjoint solution IT( to) of 
the integral Riccati equation is a nuclear operator. 

Proof Defining the bounded operator Ct 0 : X —> 1 2 (to, b: U x Y) by 

(Cto x o)(0 = ( F ^ L{ . } ) T d;to)x 0 , L = F-'B* H. 

Cf 0 is Hilbert-Schmidt by Theorem 12.31 1 The second IRE (]4]) can be rewritten as 


n(t 0 )x = Tl(b,t 0 )GT L (b,t 0 )x + C* Q C to x, x G X. 

Form Theorem l2.3l 3 and the nuclearity of G, \\(t) is a nuclear operator. □ 

Lemma 2.5. Assume T(-, •) and Tf-, •), i G N, are mild evolution operators which are 
uniformly bounded by A t, Di,D G L£°(a, b]X, X) satisfy || Di(t)x — D(t)x || —>■ 0 as 
i oo for every x € X and sup,j{||i9,: ||oo, H-ZTHoo} ^ Ad- Tjyfi-,-), Tpf,-) denote the 
perturbed evolution operators corresponding to the perturbation ofTfi-, •) by Di and T(•, •) 
by D. If || T. L (t, t)x — T(t, r)x|| — > 0 as i —> oo for x G X, then for any (t, r) G T^ and 
x G X, 

\\T Di (t 1 T)x - T D (t,r)x\\ -A 0, i —> oo. 

Proof As in O, we construct T Di (t, r) as T Di ( t, r) = Yl™=o T Di,n(t, r), where 


TDi,o{t,r) 


m,r), 



Tfit, s)Di(s)T Ditn _i(s,r)xds, x G A". 


By induction we obtain \\TD i>n {t, r)|| ^ ^ • To(f, r) can be con¬ 

structed in a similar manner with the same upper bound. 

Defining d i)n (t, r) = T Dun (t , r) - T D)U (t , r), we have r) = Tfit, r) - T(t, r), 


di, n (t,r) = J Ti(t, s)Di(s)di, n -i(s,T)ds + j Tfit, sfiDfis) - D(s)]T d ^i(s, r)ds 
+ J [Ti(t,s) - T(t, s)]D(s)T D: n _i(s,r)ds. 
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The uniform boundedness of {T/y (t. r) }, e i\ and Tp>(t, r) implies 


Y^di, n {t,T)\\ ^ ||T A (f,r)|| + \\T D (t, r)|| < oo. 

71=0 


Due to supj{11 Di 11ooi PHoo} ^ Ad and T(-, •), Tfi-, •) are uniformly bounded, the mild 
evolution operators Tb(-, •), Tyy(-, •), are uniformly bounded, and further for any n £ IN, 
supj sup,- t r ) er b || di tn (t,T)\\ < oo. Meanwhile, since \\Di(t)x—D(t)x\\ —y 0, ||di,o(^i 'r)rc|| 
|| Ti(t, t)x — T(t, r)x|| —y 0, i —>• oo. Hence, 


||dy n (f,r)a;|| ^ J \\Tfit, s)||||A(s)||||dy„_ 1 (f,T)x||ds 

+ J \\Ti(t, s)||||[A(s) - D(s)]T D , n -i(s,T)x||ds 
+ J ||(Ti(f,s) - T(t,s)){s)T Di n_i(s,T)x\\ds—> 0, i ^ 00.(5) 
By dominated convergence theorem, 

OO 

\\T Di (t,r)x - T D (t,r)x || ^ ^2 \\<k,n(t,T)x || ->• 0, i ^ oo. 


71=0 


□ 

Corollary 2.6. For any mild evolution operator T(-. •) with uniform bound Ay and Di, D £ 
Lf(tQ,b\X,X) with sup j {11 Di 11 oo > Halloo} ^ A D , if\\Di{t)-D(t)\\ -A 0, thenforT D f, •) 
which is the perturbation of T(•. •) /yy /T, mid T/y (•, •) which is the perturbation evolution 
operator ofTf, •) Ty D, we have 

\\T Di {t,T) -T D {t,r)\\ ->-0, i -A oo. 

Proof. From the assumptions, let 7) = T in Lemma 1231 replace Q by 


l|d*,n(i,r)|| ^ 


J ll T ( i > s )llllA(s)||||di,n-i(t,r)||ds 

+ [ ||T(i, s)IIHA(s) - D(s)||||Tnn(s,r)||ds —> 0, i ^ oo. 


Then, we can prove the uniform convergence of T/y (i, r) by the dominated convergence 
theorem in the similar way with Lemma [231 □ 

Theorem 2.7. We consider the time-varying system O with the location-dependent input 
operators and the cost functional ©• Assume {B r } rG fim satisfies lim r _> ro ||T r — B ro = 
0 , for some rg £ f l m , U and Y are finite-dimensional and G is a nuclear operator, then 


lim ||n r (i) — n ro 

i—^ro 


1=0, t £ [to, b] 


and there exists an optimal location r such that 

h{t) = ||IIf(io)||i = inf ||n r (t 0 )||i- 

r£il rn 

Proof. Similar to Theorem 12.21 there exists 5 > 0 such that sup reB ( r0i( 5 ) | B r \\ < oo,ro G 
Q rn and for every x G X and t G {to,b\, 

U r (t)x -► U ro (t)x, r —> ?’o- 

Further, from (0. we have II,. are uniformly bounded with An for any t G [to ? b] and 

r G B(r 0 ,<5). 

Defining the operator C t)T : X —> L 2 (T, b;U x F), t G [fo, 6], 


(Ct, r ®(i))(-) 


C(-) 

-F5(.)i?;(-)n r (.) 


T L:r (-,t)x(t). 


Corollary I2.4l has shown C/ . r is a Hilbert-Schmidt operator and 


n r (f) 


TZ ir (b,t)GT L , r (b,t)+Cl r Ct, r . 


( 6 ) 


is nuclear if G is nuclear. 

Now let us show that Ct, r uniformly converges to C //ro in Hilbert-Schmidt norm. Let 
{ e *}f=i anc i {G.}fL\ be respectively the orthogonal basis of U x Y and X, then 

II Ct,r C t,ro H/tS 

°° pb P+1 

— ^ ^ ^ (a) {Ct,ro t'i) (a) j 67 ) U xY ds 

i=l Jt j =1 

°° ,>ft P+9 

j=i 

-T£ >ro (a, t)[C*(a), L* ro ( s)F~ 2 (s)}ej} x \ 2 ds 
b oo p+q 

EE I (et,r£ !r ( S ,f)[C*( S ),L;( S )F^(s)] ej 

*=i j=i 

-r£, ro (a, t) [C" (a), L; o (a)F 3 (a)]e,>* | 2 da 
p+q 

- £ 

i=i 

where L r = 1 

From Theorem l2.21 we have lim r _j. ro ||L r (t) — L ro (t)\\ = 0 and ||L r Hoo < oo. Then, 


r£ |r (a,i)[C*(a),L;(a)F3( a )] ej .-3- iro ( s ,i)[C'*(a),L^(a)i?’3(a)] ej 


\\xds, 



lim \\B r (t)L r (t) - B r (t)L r 

r—o 


= 0 
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and H-Bj-Lj-lloo < oo. Hence, from Corollary 12.61 for any (s, t) G r£ , Tj j.{s. t) uniformly 
converges to Ti,^ 0 (s,t). Therefore, 

II - r£ >7 . 0 ( S ,t)[C'*( S ),^ 0 ( S ) J P5( S )]e J ||^ 

< ll(^L,r( s ) t) ~ T L>ro (s,t)*)[C*(s), L*(s)F^(s)]ej}\\ 

+ ll T L,ro( S >0[°:( L *( S ) - L r 0 ( S )) F ^( S )] e i]ll -► °) r r 0- ( 7 ) 

By dominated convergence theorem, || Ct, r ~ Ct,r 0 \\ns —t 0, r —> tq• Further, if G is 
a nuclear operator, 

||n r (t) - n ro (t)||i 

< ll^, r (6,i) - r£ iro (6 > t)||||Gr L , r (6,t)|| 1 + \\Tl ro (b, t)G\\i\\T Lj1 .(b, t) - T L , ro (b,t)\\ 

+ II C*,r ~ Ctr 0 \\Hs\\Ct,r\\HS + ll^ro W^S\\Ct,r ~ C^oWhS 0, r ->■ r 0 . 

By the compactness of Q m , the optimal location f exists in nuclear norm. □ 


3 Convergence of optimal control locations 

In practice, the integral Riccati equation in an infinite-dimensional space cannot be solved 
directly. Usually, we approximate and solve it in finite-dimensional space by a sequence of 
approximations from various numerical methods. Let X n be a family of finite-dimensional 
subspaces of X and P n be the corresponding orthogonal projection of X onto X n . The 
finite-dimensional spaces {X n } inherit the norm from X. For every n G IX, let T n (-, •) 
be a mild evolution operator on X n , B n (t ) G L£°(fo, b; U, X n ) and C n (t) = C(t) P„, 
G n € C(X n ). This defines a sequence of approximations 

x(t) = T n (t,t 0 )x(to) + [ T n (t,s)B n (s)u(s)ds, tG[t 0 ,b] 

■Jto 

with the cost functional 

J n (t,x,u ) = (x(b),G n x(b)) + J (C n (s)x(s),C n (s)x(s)) + (u(s),F(s)u(s))ds. 

We denote the optimal control of the approximation by u n (t) = —L n (t)P n x(t), t G 
[to, b], where L n (t) = F~ 1 B*Tl n , the perturbed evolution operator of T n (-, •) by —B n L n 
by ■) and the Riccati operator of the approximation by Fl rt . 

In order to guarantee that n n (f) converges to n(t), the following assumptions are 
needed in the approximation of control problem for partial differential equations | T51 l . For 
each x G X,u G U,y Gf, when n —> oo, 

(al) (i) T n (t, s)P n x T(t, s)x\ (ii) T*(t, s)P n x -J- T*(t, s)x 

and sup n ||T n (t, s)|| < oo, (t,s) G V b tQ . 

(a2) (i) B n (t)u —>• B(t)u\ (ii) B*(t)P n x —>• B*(t)x, a.e.. 
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(a3) (i) C n (t)P n x C(t)x ; (m) G*(f)y ->■ C*(t)y, a.e.. 

(a4) sup n ||G n || < oo and G n P n x —>■ Ga;. 

Before we study the uniform convergence from n n (f) to 11(f), we study under which 
condition the compactness of 11(f) can be guaranteed. The following lemma shows this. 

Lemma 3.1. We consider the time-varying system 0 with the cost functional 0. If B(t), 
C(t'), t £ [fo, b] and G are compact operators, then the unique solution 11(f) of the integral 
Riccati equation © is compact. 

Proof Denote S = C*C + UBF^B* II, 

n(f) = lt(b,t)GT L (b,t) + £ TZ(s,t)S(s)T L (s,t)ds. 

Since B(t), C(t) and G are compact, Tf{b,t)GTi,{b,t) and T£(s, t)S{s)Ti,{s, t), (s,t) £ 

T'f arc compact. Let us only consider the integral part of 11(f) firstly. It is clear that there ex¬ 
ists a set of orthogonal projections {P,,} to some finite-dimensional spaces X n , n G IN such 
that Hindoo \\P n Tl(s,t)S(s)T L (s,t) - T£(s,f)S'(s)T L (s,f)|| = 0. Then, since T L ( V ) 
and S(-) are uniformly bounded, it is easy to obtain P n T£,S'T^ is also uniformly bounded 
in any time and n. By the dominated convergence theorem, 

lim || f P n Tf(s,t)S(s)T L (s,t)ds — [ Tf(s, t)S(s)T L (s, f)ds|| = 0. 

n ^°° Jt Jt 

Obviously, fj 1 P n Tf(s, t)S(s)T[f s, t)ds is still finite-rank operator and bounded, so it 
is compact. 

Therefore, jj’ T£(s, t)S(s)Tf j (.s, t)ds is compact. Further, fl(f) is compact. □ 

The following theorem shows the uniform convergence of II„ (f). 

Theorem 3.2. For the sequence of approximations under the assumptions (al) — (a4), if 
B(t), C(t), t £ [fo, b] and G are compact operators and linin.^oo \\B n — Pn-BHco = 0, then 

lim ||n n (f)P n - n(f) || =0, f € [f 0 , b\. 

n—>oo 

Proof From lim^oo \\B n - P n .B||oo = 0 and sup ie[to fc ] ||5(f)|| < oo, we have 

sup ||S„(f)|| < oo. 

nElN ,£E[to M 

Moreover, because B(t) is compact and P n is strongly convergent to the identity operator 
I, lim, woo \\P n B(t) - B(t)\\ = 0, f £ [f 0 , b\. Further, 

\\B n (t) - B(t)\\ ^ ||5 n (f) - P„5(f)|| + ||P n S(f) - B(t)\\ —>• 0, f £ [f 0 ,6], n —> oo. 

Meanwhile, by the uniform boundedness of ||T n (-, -)||, ||G n ||oo and ||G n || and [ T5l The¬ 
orem 5.1], for any x £ X, 

lim ||II n (f)x — n(f)s|| =0, f £ [f o, b] - 

n—>oo 
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Similar to the proof of the uniform boundedness of II r in Theorem 12.21 for the ap¬ 
proximations with arbitrary feedback control u n (t) = L n (t)x(t) = L(t)P n x(t), L G 
L^°(fo, b-,X,U), there exists An > 0, such that sup n Hllnlloo < An- 

To proof the uniform convergence of II n (t), we define S n = C*C n + \l n B n F ^ B*JA„ 
and S with the similar way, then 

||n n (f)P n -n(f)|| ^ \\(Tl n (b,t) — Tl(b,t))G n P n \\\\T Ln (b,t)P n \\ 

+ \\T* L (b,t)\\\\(G n Pn - G)T Ln (b,t)P n \\ + \\TZ(b,t)\\\\G(T Ln (b,t)P n - T L (b,t))\\ 

+ J II [T* Ln C a,t ) - T* L (s, t)}S n (s)P n || ||T in ( S , t)P n \\ds 
+/r£( S ,f)iiii5 n ( S )p n -5( S )iniT Ln ( S ,f)p n ||d S 
+ J* II T L (s, t)|| II S(s)(T Ln (s, t)P n - T L (s, t))\\ds. 

As aresult of the uniform boundedness of ||T n (-, -)ll> l|n ?l ||oo and ||f?n||oo i n n, IIZvjt. ||oo 
is uniform bounded and 

||L*(t)-P n L*(f)|| < ||F- 1 || 0O (||n n (f)||||S n (f)-P n S(f)|| 

+ ||n n (f) - P n II(f)||||B(f)|| —> 0, r —>• r 0 . 

so Hindoo \\L n (t)P n - L(t)\\ = 0 and linr^^ \\B n (t)L n (t)P n - B(t)L(t)\\ = 0. 
According to Lemma 1231 and assumption (al), 

lim || T Ln (t, s)P n x - T L (t, s)x|| = 0, 

n—too 

lim \\Tl (t,s)P n x-Tl(t,s)x\\=0, x € X. 

n—too 

Finally, because of the compactness of the self-adjoint operator G n and G, we have 

lim || (T Ln (t,s) - T L (t,s))G n P n \\ = 0 

n—too 

and lim^oo \\G(T Ln (t,s)P n - T L (f,s))|| = 0. Meanwhile, H^Hoo < ||C'*C' ri ,||oo + 
||n n 5 n F- 1 s*n n || 00 < oo, n G IN. Since C n = CP n is compact, 

\\S n {t)P n -S(t)\\ 

^ II C* n (t)C n (t) - C*(t)C(t )II + \\L* n (t)F(t)L n (t)P n - L(t)F(t)L(t )|| 

^ ||C*(f)P n - C'*(t)||||(7 n || 0O + ||C , *|| 00 ||C'(f)P n - C(t)\\ 

+ || L* n (t) - L*(t)||||F||oo||£n||oo + ||L*|| 00 ||F|| 00 ||L n (f)P re - L{t )|| -»■ 0, n -A oo. 

By dominated convergence theorem, Ulln^fjPjj — II(t)|| —> 0, n —> oo. □ 

Next we show that the optimal control locations of approximations converge to the 
optimal control location of the original system. 
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Theorem 3.3. Under the assumptions (al) — (a4) and further assume B rn = P n B r , r G 
Ll m , if B r (t), C{t ) and G, t G [to, 6] are compact operators and lim,.-^ \\B r — B ro || = 0, 
then 

i n (t)^ri(t), f n —y r, n ->• oo. 

Proof From Theorem 13.21 lim n _ > . 00 [|n r>n (t)P n — n r (t)[| = 0 , r G Ui m . 

Since lim r .—||FA> -I^olloo — 0, 

||-®r,n -®ro,n||oo ^ ||Pra 1111 B r -B^qUoq f 0, V > Tq. 

From Theorem 12.21 for any n G I, there exists l n (t) = inf re Qm ||n rin (i)||. 

On one hand, 

(■nit) = inf ||II rn (f)|| ^ ||IIf,„(t)|| < ||n.r, n (f) — Ilr(t)|| + 11Ilf.(t)|j 

rEi2 m 

—> 11 Ilf (t)\\ = l(t), n —y oo, 

so limn^oo sup n l n (t) ^l. 

On the other hand, there exists a subsequence {l, lk (t) \ such that lini/._ >oc i nk (t) = 
lim^-^ooinf n £ n (t), where i nk (t ) = inf rG n™ |[II r)nfc (t)|| = ||n rnfcinfc (t)||. Due to the com¬ 
pactness of Q"\ without loss of the generality, we assume lim/,^^ f nk = r, 

II ^^n k ,n k I loo ^ ||Prifc 111| ||oo T ||Pnt, Bf Hf ||oo t 0; k y OO 


and 


\\Bfn k ,n k it) n r (f)|| ^ l|Hf nfc ,n fc (i) n^(f)|| + ||n r „ fc (f) n r (t)|| to, k—y oo. 

( 8 ) 

Hence, 

lim inf £ n (t) = lim £ nk (t) = lirn ||IIf „ fc (i)|| = ||n f (t)|| > ||n.f(f)|| = l r (t), 

rwoo n fc->oo k^too k 

so lim n ^. 00 i' n (f) = l(t). Further, lim n ^. 00 i n (f) = lim^ooinf n £ n (t) = £(t), so 

lim Pf nk (f)|| = ||n f (f)|| = ||n f .(t)||. 

k—yoo K 

By the continuity with respect to r nk in ©, lim*,-^ f Uk = r. □ 

For the proof of the uniform convergence of the Riccati operators of the approximations 
in nuclear norm for stochastic systems, we need the following lemma. 

Lemma 3.4. Let G be a nonnegative nuclear operator in a separable Hilbert space X and 
assume that T n strongly converges to T, T n , T G C(X) are uniformly bounded by A t- Then 


lim \\(T n — T)G||i = 0. 

n—>-oo 
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Proof. Assume {e,} is the orthogonal basis in X and there exist a partial isometry V such 
that G = V\G\, where |G| = (G*G)i, then, 

||(T n - T)V\G\hi\\ < \\T n - T\\\\V\G\hi\\ ^ 2A T ||P|G|^||. 

Because of the strong convergence of T n , lim^^oo || (T n — T)V|G| 2 e* |j = 0. 

Since G is a nuclear operator, then |G| 2 is a Hilbert-Schmidt operator, so 


£ ll( r - - T ) v \ G \ he i\\ = 2AtE II^GI^I 


< 00 . 


i=1 


i =1 


By the dominated convergence theorem, 


lim \\(T n -T)V\G\2\\ HS 

n—too 


lim V ||(T n - T)V\G\hi\\ = V lim ||(T n - T)V\G\hi 

r?,-4oo zJ zJ n.—±on 


= 0. 


i— 1 


i=l 


Then, ||(T n — T)G||i < ||(T n - T)V\G\ 5 ||^ 5 || |G| 2 \\ HS 0, n 00 . 


□ 


Associated with Corollary 12.41 the following theorem guarantees the uniform conver¬ 
gence of the Riccati operators of approximations to the Riccati operator of the original 
system in nuclear norm. 

Theorem 3.5. For the sequence of approximations under the assumptions (al) — (a4), if 
U and Y are finite dimensional, lim n _ > , 00 \\B n — P n -B||oo = 0, G is nuclear operator and 
lim^oo ||G n P n - GHi = 0, then 


lim ||II n (t)P n - n(t)||i = 0. 

n—too 

Proof Defining Ct in the same way with Corollary I2.4l and define C Ln by substituting n into 
r in ©, from Theorem l2.3l l. II n (t) = Tfi, (6, t)G n Tr, n (b, t) +C* n C Ln is nuclear. The same 
with Theorem 13.21 we also have the uniform boundedness of ||T n (-, -)||, ||n n ||oo 11 Ti11oc, 
HL^Hooinnandlim^oo \\L n (t)P n - L(t)\\ = 0, lim^oo \\T Ljl (t,s)P n x-T L (t,s)x\\ = 
0. Hence, similar to Theorem 12.71 

+ H ? zXM)[C' n (s) - G(s), (L*(s) - L*(s))F^(s)]ej]\\ds —* 0, r -> r 0 , ( s,t ) € T b to . 
Then, since G is nuclear operator with lim,,^,*, ||G n P n — G|| 1 = 0, 

||n n (f)p n -n(f )|| 1 

< \\TUb,m\\Gn(T Ln (b,t)Pn - T L (b,t)) 111 + II(T£ n (M) - T£(M))G n ||i||T L (M)|| 
+ \\Tt(b,t)\\\\G n Vn - G||i||T L (6,t)|| + ||C* n - C* t \\ H s\\Ct,n\\ H S 
+ ||G t *||ifs||Ct )n — C t \\HS —> 0 , n —> 00. 


□ 
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Theorem 3.6. Under the assumptions (al) — (a4) and further assume B rn = P;, B r , 
r £ X m , if the input space U and the output space Y are finite dimensional, lim r _> ro \\B r — 
B ro | = 0, G is nuclear operator and lim n _ > . 00 || G n P n — G||i = 0, then 

h,n(t) ->• £i(t), r n ->r, n —>■ oo. 

Proof From Theorem 12.71 and Theorem 13.51 we have lim r _ >ro ||n r (t) — n ro (t)||i = 0 
and lim^^oo ||H, v , t (t)P n — II r (i)||i = 0. The same with Theorem 13.31 we have £i >n (t) ^ 
||IIf (t)||i = £\(t), n —> oo. Besides, there exists a subsequence {in ,n k (t)} such that 

lim inf l\ n {f) = lim h,n k (t) = lim ||II f n (f)||i = ||n f (f)||i ^ ||II^(i)||i =h(t). 

n^roo n fc—roo oo 

Therefore, lim^oo h,n(t) = h (*) and lim*..^ ||n^ fc)7lfc (i)|| 1 = ||n ? (i)|| 1 . By the 
continuity in Theorem l2.71 lim^oo r nk = r. □ 

4 Kalman filter in Hilbert spaces and the duality of LQ optimal 
control problem 

There are several literatures 0, 0, lll3l . IT8l discussing the Kalman filter in different 
approaches and the duality to the linear-quadratic optimal control. Flowever, to author’s 
knowledge, these derivations involve the generator of semigroups or evolution operators. In 
this section, without the differentiability of evolution operators, we first derive the Kalman 
filter in real separable Hilbert spaces. Further, we will discuss the duality between Kalman 
filter and linear-quadratic optimal control. 

Let (ft, B. p) be a complete probability space and X, £,y be a real separable Hilbert 
spaces. First, we define some basic concepts of probability theory in Hilbert spaces 0, 

m. 

Definition 4.1. The map x : ft -X X is a X— valued random variable if it is strong measur¬ 
able with respect to a measure //. 

Definition 4.2. // is a totally finite measure on X if for any X—valued random variable x, 
J (i ||x||d/r < oo. Further, if there exists x £ X such that 

(x,h) = E (x,h) = / (x,h)p(dx), Vh £ X, 

Jn 

x is called the mean or expectation of x and denoted by Ex. 

Definition 4.3. For any X—valued random variable x with mean Ex, the covariance oper¬ 
ator P of x, also denoted by Cov(x), if it exists, is given by 

(Phi, hf) = (hi, Phf) = / (x - Ex, hi)(x - Ex, h 2 )p(dx), \Jh\,h- 2 ^X. 

Jn 

Definition 4.4. The random variables x, y whose expectations exist are independent if 
E«s, 2/» = (E(x),E(y)). 
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Definition 4.5. Let fibe a probability measure on X. If for any x £ X, the random variable 
(x, ■) has a Gaussian distribution, then // is called a Gaussian measure. Further, we denote 
x of the Gaussian measure with mean x and covariance F by x ~ N(x, P). 

Definition 4.6. t £ IR} is a set of white noises if for any t £ [0, +oo], there exists a 

covariance operator W (t) such that ui(t) ~ N{ 0, W (t)) and for any t s, uft) and uj(s) 
are independent. 

We consider time-varying systems on Hilbert spaces given by 

x(t) = M{t,t 0 )x(t 0 ) + f M(t, s)[B(s)u(s) + D(s)u(s)]ds, (t,t 0 )eT^ o , (9) 

Jt 0 

where M(-, •) is a mild evolution operator on X. x(t) and u(t) are random variables with 
values in X and £, respectively and oj{t) ~ W(0, W(t)) is the white noise. Further, we as¬ 
sume u £ L 2 (t 0 , b ; U ), B £ Lf (t 0 ,6; U, X), B* £ Lf (t 0 , 6; X, U), D £ Lf{t 0 ,6; £, X). 
We consider the following observation system 

y(t) = H(t)x(t ) + E(t)u(t), t £ [to, 6], (10) 

where H £ L“(fo, b ; X, y), E £ o, b ; £, y), y(t ) and u{t) are random variables with 
values in y and £, respectively and v{£) ~ N(0,V(t)) is the white noise and V(t) is a 
coercive operator.. 

In our paper, we only consider the integral form of time-varying systems. Let Y t = 
(y(s), to ^ $ f t}. the linear unbiased estimation of the filter problem x(t\t) of x(t) lfl9l 
has the form 

x(t\t) = M(t,to)x(t 0 \t-i) (11) 

+ / M(t, s)B(s)u(s)ds + / Kf(t,s)[y(s) — H{s)x(s\s)]ds, 

J to Jto 

where x(to\t-i) = E(x(to)), Cov{x{tf)) = P(to\t-i) and Kff, ■) is an unknown linear 
gain operator. 

Denoting x(t\t) := x(t) — x(t\t), P(t\t ) := Cov(x(t )) and R(t) := E(t)V(t)E*(t), 

we obtain the following theorem 

Theorem 4.7. For the time-varying system © with the observation system CO]), the linear 
unbiased estimation of the filter problem x{t\t) ofx{£) is optimal if the linear gain operator 
in (fill) is given by Kf(t,r) = M(t, r)P(r|r)iF*(r)i? _1 (r), t f t. 

Proof By Wiener-Hopf’s equation Ifl3l . | |T9| |, xit\t) minimizes the minimal covariance if 
and only if E(£(t),/ii)(y(r) — H(t)x(t\t), hf) = 0, r < t, h\, hi £ X. Further, 
according to ||6l Corollary 6.3], E(x(t\t), hi)(x(t\t), hf) = 0. Hence, on one hand, 

E(x(t), hi){y{r) - H(t)x{t\t), h 2 ) 

= E (x(t), hi)(H( t)x(t\t), h 2 ) + E (x(t), hi)(E(T)i/(T), h 2 ) 

= E(M(t,r)x(r),/ii)(iF(r)F(r|r), h 2 ) - E(M(t,r)x(r|r),/ii)(Ff(r)x(r|r), h 2 ) 

= E(M(f, r)x(f|t), hf){H{ t)x{t\t), h 2 ) = {hi, M{t,T)P(r\T)H*(r)h 2 ). 
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On the other hand, 

E(x(t|t), hi)(y(r) - H(t)x(t\t), h 2 ) 

= E( [ K f (t,s)[H(s)x(s\s) + E(s)u(s)]ds,hi)(y(T) - H(t)x(t\t), h 2 ) 


= E( 
= E( 


'j{ j Kf(t, s)E(s)v(s)ds, h\)(H (t)x(t\t) + E(t)u{t), h 2 ) 

)(J Kf(t,s)E(s)u(s)ds,hi)(E(T)i/(T),h 2 ) = {hi, K f (t, t)R(t)}i 2 ) 

Therefore, Kf(t,r)R{r ) = M(t, t)P(t\t)H*(t). Since R(t) is coercive, we obtain 
Kf(t,T ) = M(t,T)P(T\T)H*(r)R~ 1 (T), t < t. 

If t = t, by the strong continuity of Kf(t, •), Kf(t,t) = P(t\t)E[*(t)R~ 1 (t). 

Defining K(t ) := Kf(t,t ) = P(t |i)i7*(f)i? _1 (i), Theorem 14.7 1 implies that 


'to 


t rt 



Mx(t, r])K{rj)H(rj)M(jj , s)K(s)H(s)x(s\s)drjds 
+ f Mx(t, s)[D(s)u(s) — K(s)E(s)u(s)]ds 

J to 

( f MK(t,ri)K(r])H(r})M(r),s) (D(s)lu(s) — K(s)E(s)u(s)) drjds 
J tn J s 


'to Js 
rt 


□ 


x(t\t) = M(t,to)x(to\t-i) — / M(t,s)K(s)H(s)x(s\s)ds 

Jt 0 

+ f M(t, s)[D(s)lj(s) — K(s)E(s)v(s)]ds. (12) 

Jto 

Theorem 4.8. Equation (1121) is equivalent to 

x(t\t) = M K (t,t 0 )x(to\t-i) + f M K {t,s) (D(s)uj(s) - K(s)E(s)i/(s)) ds, (13) 

Jt 0 

where Mx{t,r)x = M(t,r)x — f* Mx(t,s)K(s)H(s)M(s,r)xds, (t,r) € T^ 

Proof. From (fl2l) . 

x(t\t) 

= M K (t,t 0 )x(t 0 \t-i) + [ M K (t,s)K(s)H(s)M(s,t 0 )x(t 0 \t- 1 )ds 

Jt 0 

— / Mx(t, s)K(s)H(s)x(s\s)ds 
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□ 


For finite-dimensional systems, the trace of the covariance of x(t\t) is considered as an 
evaluation of the estimation errors. For systems on Hilbert spaces, similarly we consider 
the nuclear norm of the covariance of x(t\t). Defining Q(t) := D(t)W(t)D*(t), we obtain 
the following theorem. 

Theorem 4.9. The covariance (if exists) of x(t\t) satisfies the IRE 


P(t\t) = M K (t,t 0 )P(to\t-i)M* K (t,t 0 ) 


(14) 


+ [ M K (t, s ) [Q(s) + P(s\s)H*{s)R- 1 (s)H(s)P{s\s)} M* K (t , s)ds. 

Jto 

Proof For x(t\t) in (IT3l) . assume its covariance P(t\t) exists and define Qt- L 2 (to,t;£ x 
S) -)• X by 


Qt 



ds. 


Its adjoint operator Q*: X —> L 2 (to, t:£ x £) is given by 



Then, we obtain 


E(x(t\t),h 1 )(x(t\t),h 2 ) 

E(M K (t,t 0 )x(to\t-i),hi){M K (t,t 0 )x(t 0 \t-i), h 2 ) 



{M K (t,t 0 )P(t 0 \t-i)M^(t,t 0 )h 1 ,h 2 ) + ( Qt.Cov 
(M K (t,t 0 )P(t 0 \t-i)M^(t,t 0 )h 1 ,h 2 ) 
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+ ( [ M K {t, s ) [Q(a) + K(s)R(s)K*(s)} M* K {t , s)Ms, h 2 ), 

J to 

Hence, for any x € X 


P{t\t)x = M K (t,to)P{to\t-i)M^(t,to)x 

+ [ M K (t, s ) [Q(s) + A(s)A(s)A*(s)] Mx(t, s)xds 
Jto 

= M K (t,t 0 )P(to\t-i)M^(t,t 0 )x 

+ [ M K (t,s ) [Q(s) + A(s|s)A*(s)A~ 1 (s)A(s)P(s|s)] M* K {t,s)xds. 

Jt 0 

□ 


A comparison with the main results of the linear-quadratic optimal control problem in 
Section [2] yields: By observing the similarity between (fldl) and the second integral Riccati 
equation related to the linear quadratic optimal control problem, it is clear that to consider 
the covariance of x(t |f) of the time-varying system ([9]) with the observations (fTOl ) is equiv¬ 
alent to consider the Riccati operator 11(6 — i) in © corresponding to the time-varying 
system 

x(t) = T(t, to)x(to) + / T(t,s)B(s)u(s)ds. 

Jt 0 

with the cost functional 


J(t, x, u) = (x(b), Gx{b)) + J (C(s)x(s), C(s)x(s)) + (u(s), F(s)u(s))ds, 

where T(t, s ) = M*(b — s,b — t), B(s ) = H*(b — s), G = P(io|i-i), C(s) = Q^(b — s), 
F(s) = R(b-s),(t,s)e V b tQ . ' 

Then, by the duality between the linear quadratic control problem and Kalman filters, 
Corollary 12.41 implies the following condition to guarantee the existence and nuclearity of 
P(t\t). 

Theorem 4.10. For the time-varying system © with the observation system d, if £ and 
y are finite dimensional and P(fo|f-i) is a nuclear operator, then the covariance of x(t\t) 
based on Yj satisfying cs exists and is a nuclear operator. 


5 Kalman smoother in Hilbert spaces 

In this section, we study the optimal linear unbiased estimation of x(t) based on Y t by 
x(r|6), t f t. We still constrain the linear estimation of x{r\t) has the form 

x(r\t) = j K s (t, s)[y(s) - H(s)x(s\s)]ds, r^t, (15) 

Jt 0 
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where K s (■. ■) is an unknown linear operator. 

Since in the case r = t, (fl5T ) with the minimal covariance is equivalent to the opti¬ 
mal linear unbiased estimation based on Kalman filter, in order to determine the optimal 
estimation of x{r\t), r ^ t, we can rewrite (fl~5l) as 

x(r\t) = x(t\t) + J K s (t, s)[y{s) — H{s)x(s\s)]ds. (16) 

Theorem 5.1. For the time-varying system © with the observation system COl), the linear 
unbiased estimation of the filter problem x(r\t) ofx{r) is optimal if K s {-, •) in dill) is given 
by 

K s ( t, rj) = P(t\t)Mk(t),t)F[*{ p)R~ l {p), r < p ^ t. 

Proof By Wiener-Hopf’s equation fl9l . lfl3l . E{x{r\t), hi){y{p) — H{p)x{p\p), h 2 ) = 
0, h\ € X, h 2 € y, for any p < t. It is clear for any p < t, E(x{r\t) 1 hi){y{p) - 
H(p)x(p\p), h 2 ) = 0 holds. Now we assume r ^ p < t. On one hand, 

E(x(r), hi)(y(p) - H{p)x{p\p), h 2 ) 

= E (x(t) - x(r\p), h 1 ){H(p)x{p\p),h 2 ) 

= E(x(t\t) - J I< s (T,s)[y(s ) - H(s)x(s\s)}ds,hi){H(p)x{p\p),h 2 ) 

= Ei(x(t\t), h\)(H( p)Mk{ p, t)H*(p)x(t\t), h 2 ) = {hi, P{t\t)M^{p, t)H*( p)h 2 ). 

On the other hand, 

E(x(r|i), hi)(y(p) - H{p)x{p\p), h 2 ) 

= E{x(r\p), hi)(H{p)x{p\p) + E{p)v{p), h 2 ) 

+E( [ K s (r,s)[y(s) - H(s)x(s)]ds,hi){y(p) - H{p)x{p\p), h 2 ) 

Jri 

{J K s {s,T)E{s)u{s)ds,hi){E{p)u{p),h 2 ) = {hi, K s {t, p)R{p)h 2 ). 


= E 


By the coercivity of R{t), we obtain K s (r,p ) = P{t\t)M^{p,t)H*( p)R 1 (p)- 
Defining x{r\t ) = x(t) — x{r\t), Theorem I5T1 implies 


□ 


x{r\t) = x{t\t) — P(t\t) J AI^{s,t)H*{s)R i {s)[y{s) — H{s)x{s\s)]ds. 

Thus, its covariance can be derived by 
Theorem 5.2. The covariance (if exists) ofx{r\t), {t, r) e r ?o is 
P{r\t)x = P{t\t)x 

— P{t\t) J Mk(s, t)H*{s)R~ 1 {s)H{s)Mk{s, T)P{T\r)xds, x G X. 


(17) 


20 


Proof. Denoting the covariance of x(r|f) by P(r|t), we obtain 


(hi,P(T\t)h 2 ) = E(x(r|f),/ii)(x(r|f),/i 2 ) 

= E(x(t|t) - P{t\t) J Mj < (s,T)H*(s)R^ 1 (s)[y(s ) - H(s)x(s)]ds, hi)(x{r\t), h 2 ) 

= E(x(t|t),/ii)(x(t|t) 

-P(t|t) j M* K (s,T)H*(s)R~ l (s)[y(s) - H(s)x(s\s)]ds, h 2 ) 

= E(x(t|t),/ii)(x(t|t),/i 2 ) 

— E(x(t|t),/ii)(P(t|t) J M^(s,T)H*(s)R~ 1 (s)H(s)x(s\s)ds, h 2 ) 

= (/ii,P(r|r)/t 2 ) - F(x(t|t),/ii)(P(t|t) J Mj < (s,t)H*(s)R~ 1 (s)H(s) 

■ (m k (s,t)x(t\t) + J M k (s, rj) (. D(ri)u(ji ) - K(ji)E(rj)o(ri )) ds,h 2 ) 

= (/ii,P(r|r)/t 2 ) 

— E(x(t|t), /ii)(P(r|r) J Mk(s,t)H*(s)R~ 1 (s)H(s)M k (s, r)x(r|r)(is, /t 2 ) 

= (/ii,P(r|r)/t 2 ) - (/ii,P(r|r) J M^(s, t)H*(s)R~ 1 (s)H(s)M k (s, T)P{r\r)h 2 ds). 
Hence, for any x € Af, (i, r) G rf Q , we get 

P(r|i)x = P(t|t)x — P(r|r) f Mffs, t)H*(s)R~ 1 (s)H(s)Mk(s , r)P(r|r)x(is. 


□ 


Theorem 5.3. For f/ie time-varying system © with the observation system CO]), if £ and 
y are finite dimensional and P(to\t-\) is a nuclear operator, then P(r\t), ( t,r) G F^. 
satisfying C3t exists and is a nuclear operator. 

Proof. By Theorem H.lOl and the uniform boundedness of Mk, H and R 1 in [to, b ], 

ll- p M*)lli < ll- p (^k)lli + Il-P(r|r)||? J ||M^(s,r)|| 2 ||P" 1 (s)||||P(s)|| 2 ds < oo, 

so P(r|i) is a nuclear operator for any (t, r) G F^. □ 


6 Optimal locations of observations based on Kalman filter and 
smoother 

In this section, we also take the observation location problem into account. The location 
parameter r is defined as in Section [2J The following theorems show the continuity of 
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P r (t\t) and P r (r|i), (t, t) E r £ in nuclear norm. For the filter problem, due to the duality 
and Theorem l2.71 we obtain the following theorem. 

Theorem 6.1. Consider the filter problem of the time-varying system © with location- 
dependent output operators and the observation system m. If H r is of the property that 
lim r ^. ro \\H r —H r o||oo = 0, £ and y are finite-dimensional, and P(to|f-i) is nuclear, then 


lim ||P r (t|t) 

r—>ro 


^r 0 (*l*)l|l =0, t£[t 0 ,b], 


and there exists an optimal location f* such that, 

£{(t) = l|-Pf/(^)lli = inf \\Pr (f|f)||l■ 

Theorem 6.2. Consider the smoother problem of the time-varying system © with the 
location-dependent output operators and the observation system ©. H r has the prop¬ 
erty that lim r .^ ro \\H r — H ro ||oo = 0. If £ and y are finite-dimensional, and P(fo|i-i) is 
nuclear, then, 

lim ||P r (r|f) - P ro (T|i)lli = 0, (^ r ) e r \ 

r—»ro u 

and there exists an optimal location r s such that for any initial time r£ [to, b], r ^ t, 
h( T \t) = ||Pp(r|f)||i = inf ||P r (r|f)||i. 

r&l m 

Proof From Lemma [5731 P r (r\t), r E Q m are nuclear operators. Flence, 

\\P r (r\t) - P ro (r|f)||i < \\Pv{t\t) - P ro (r|r)||i 

+ [ \\P ro (T\T)M^ ro (s, t)H* q (s) — P r (T\r)M^ r (s, t)H*(s)\\ 

Jt 0 

■ \\R~ 1 (s)H ro (s)M K}ro (s,T)P ro (T\T)\\ids+ f \\P r (T\T)Mj <r (s,r)H*(s)R~ 1 (s)\\ 1 

Jt 0 

■ \\H ro (s)M K . ro (s, t)P T0 {t\t ) - H r (s)M Ktr (s, r)P r (r|r)||ds, 

Since P r (t),r E £l m are nuclear operators and P _1 (f), H r (t), M^ ro (f, r) arc uni¬ 
formly bounded for (t,r) E r£ , then ||P _1 (s)P ro (s)Afx,r 0 ( s >' r )-f , r 0 (' r l r )l|i < 00 an( f so 
is its adjoint. 

By Theorem [6T] and dominated convergence theorem, we obtian 

||P r (r|f) - P ro (r|f)||i ->• 0, r —>• r 0 . 

Because of the compactness of £l m , there exists the optimal location of observations such 
that l\{r\f) = ||P^(r|f)||i = inf rG nm ||P r (r|f)||i. □ 

Next we consider a sequence of approximations of time-varying systems in order to 
study the convergence of optimal observation locations based on Kalman filter and smoother. 
Let X n be a family of finite-dimensional subspaces of X and P„ be the corresponding 
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orthogonal projection of X onto X n . The finite spaces { X n } inherit the norm from X. 

For n G IN, let M n (-, •) be a mild evolution operator on X n , D n (t) = P n D[t) and 

H n (t ) = H(t) P n , t G [fo)&]- I n order to guarantee that P n (t\t) converges to P(t\t), the 

following assumptions arc needed in the approximation of observation problems for partial 

differential equations. For each x £ X, l o € £, y € y 

(Al) (i) M n (t, s)P n x —> M(t, s)x\ ( ii ) M*(t, s)P n x —> M*(t, s)x 

and sup n || M n (t, s)|| < oo, for any (t, s ) G T b Q . 

(A2) (i) D n (t)u —> D(t)uj\ (ii) D*(f) P n x —> D*(t)x, a.e. t G [to,b], 

(A3) (i) H n (t)P n x-*■ H(t)x; (ii) H*(t)y —> H*(t)y, a.e. te[t 0 ,b}. 

(A4) P n (t 0 \t-x)P n x P(t 0 \t-i)x and sup n ||P„(i 0 |i-i)|| < oo. 

The next theorem shows the uniform convergence of the approximations of covariances 
of the Kalman filter and smoother in nuclear norm. 

Theorem 6.3. Assume that the assumptions (Al) — (.44) are satisfied. If £ and y are 
finite-dimensional, linin^oo ||P n (fo|f-i)P,i — P(to[i-i)||i = 0 and P(io[i-i) is nuclear, 
then 

lim \\P n (t\t)P n - P(t\t)\\i = 0, 

n—too 

lim \\P n (T\t)P n - P(t|*)||i = 0, (t,r) G T b t 

n—too u 

Proof Due to the duality between Kalman filter and LQ optimal control problem, according 
to Theorem [33J we have 

lim \\P n (t\t)P n - P(t|t)||i = 0, (t,r) G T b t0 . (18) 

n—>oo u 

Then, 

\\P n (r\t)P n - P(r|t)||i ^ ||P n (r|r)P n - P(r|r)||i 
+ [ ||P(r|r)M^(5,r)^( S )-P n P n (r|r)M^( S ,T)^( S )|| 

■\\R- 1 (s)H(s)M K (s,T)P(r\T)\\ 1 ds + j' \\P n P n (T\r)M^ n (s,T)H* n (s)R- 1 (s)\\ 1 
■ \\H(s)M k (s,t)P(t\t) - H n (s)M Kt n(s,T)P n (T\T)P n \\ds. 

where, according to Lemma [231 and (fT8l) . 

\\H(s)M K (s,t)P(t\t) - H n (s)MK,n(s,T)P h (t\t )P n || 

^ \\H(s)M k (s,t)\\\\P(t\t) - Pu(t\t )P n || + I! H(s) - H n (s)\\\\M K (s,r)P n (T\r)P n 
+ \\H n (s)\\\\(M K (s,T) - M K: n(s,T))P n (T\T)P n \\ 0, Tl > OO. 

and so is its adjoint operator. 

By the uniform boundedness of P(t\t), Mx(t,s), H n (t) for f G [to- Ii, we have 
||P n (r|f)P n - P(r|t)||i -A 0, n -5- oo. □ 
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Now let us take the location of observations into account and show the convergence of 
optimal observation locations of approximated covariance of Kalman filter and smoother. 


Theorem 6.4. Assume that the assumptions (A\) — (.44) are held and that H rn = H r P n 
and lim r ^ ro \\H r — H ro || oo = 0- If £ an< I y are finite-dimensional, P(fo|^-i) is nuclear 
and li mn^oo ||_P n (t 0 |i-i)Pra - P(io|i-i)||i = 0, then 


i(*)> i{(T\t), r s n ^r s , 

Proof Follows by duality and Theorem 13.61 


(t,r) £ T b to , n-> oo. 

□ 


7 Application 


As a popular data assimilation method, the ensemble Kalman filter and smoother are widely 
applied in meteorology. Flence, we consider a linear advection-diffusion model with Q := 
(0, 5) x (0,5) x (0,1) on a fixed time interval [0,3] based on the Kalman filter and smoother, 
the theoretical foundation of the ensemble Kalman filter and smoother, as an example: 


d5c ddc dSc d . d5c 

-W = ~ v ‘te- v >-at + & (A (2)+ & 

5c(t. 0 ) = 5c 0 , Se(t 0 ) = 8e 0 , 8d(t 0 ) = Sd 0 , 


Sd, 


where Sc, Se and Sd are the perturbations of the concentration, the emission rate and de¬ 
position rate of a species, respectively. v x and v y are constants and K(z) is a continuous 
differentiable function of z. 

a a 9 d 

Defining A x := -v x jfe, A y := - v y j | and D z := — (K(z) — ) with domains 

V{A X ) = {/ € L 2 (Q) | A x f G L 2 (n), /(0, y, z) = /(5, y, z)}, 

V(A y ) = {/ G L 2 (Q) | A y f G L 2 (D),/(x,0,z) = f(x,5,z)}, 

V(D Z ) = {/ G L 2 (Q) | D z f G T 2 (fl), f z (x, y, 0) = f z (x,y, 1) = 0} 


and denote by S x , S y and S z the semigroups generated by A x , A y and D z . S is the semi¬ 
group generated by A x + A y + D z with the domain V = V(A X ) n V(A y ) n V(D Z ). 

In particular, in order to include the emission rate into the state vector as optimized 
parameter, the dynamic model for the emission rate with constant emission factors llT2l is 
established as 

5eft) = M e (t, s)Se(s), 

where M e (t,s) = G L(L 2 (fil)), e^f) G L 2 (Q) is termed as the background knowl¬ 
edge of the emission rate, which is continuous in time and sup (t s ) g rg llf^y || < oo. Ac¬ 
cording to Definition 12.11 it is easy to show that M e f, ■) is a self-adjoint mild evolution 
operator. 

Ignoring the model error, the model extended with emission rate is given by 

/ 6c(t + At) \ 
y Se(t + At) J 
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M(t + At,t) ^ ) _ ( S ( t + / ^ s ^ Sd ^ ds ^ (19) 


where 


M(f , A tt)=( S ( At ) It t+At W + At ~ s ) M e (s, *)<*« \ 
1 V 0 M e {t + At,t) ) 


also satisfies Definition ^. II 

The numerical solution is based on the symmetric operator splitting technique 0, If27l 
with space discretization via finite difference method with discretized intervals Ax, Ay 
and A z in three dimensions. We assume that the grid points {r,}” =1 have the coordinates 
{ (x Ti , y ri , z Ti )} and define the projection P„ : A 2 (S >) —>■ R" 


(P nf)i ■ = 


i Jft 


f(oj)duj, i = 1 , • • • , n. 
Aj „ I Aai 


( 20 ) 


where D* = [x ri - ^,x ri + ^r] x [y ri - ^-,y Ti + x [z r< - + ^], Vj 

is the volume of H*. Defining 5 n (At) := 5 Xjn (^)S' y , n (^)5^ in (At)52 /in (4 i )5 , ;Ei n(^), 

according to |QQ Theorem 3.17], we obtain 


lim \\(S n (At))^P n f-P n S(t)f\\ = 0, f € L 2 (D). 

n—>-oo,A£—>-0 


( 21 ) 


With the same space discretization for 5c, the approximation of the emission rate is given 
by P n 5e(t) = M e ^ n (t, s)P n <5e(s), where M e ^ n {t, s ) is a diagonal matrix with the diagonal 
given by 

,. (1 . r , ,/oi e b(t,u)du> f n e b (t,uj)du 

diag (M e , n (t,s)) = (- ' 


L e b (s, u)du ’ ’ L e b {s, uj)du ' 


Then, we can easily get 

f n . e b (t, ui)du 


f n e b {s,uj)duj 


{P n f)i-(P n M e (t,8)f)i 


0, n —>• oo, / £ L 2 (D), 


so is the adjoint of M e (t, s). The extended model with operator splitting discretized in space 
can be written as 


5c n (t + At) 
6e n (t + At) 


= M n (t + At,t) 


5c n {t) 

5e n {t) 

S x ,n(^r)Sy,n{^r) ft +At S Zj n(t + At - s)5d n (s)ds 

0 


where 5c n (t) = P n 5c(t), 5e n (t) = P n 5e(t) and 5d n {t) = P n 5d{t) 

M (t + At t) = ( S x , n (^)S y , n (^)i t t+At S z , n (t + At-s)M e , n (s,t)ds \ 

n[ ’ V 0 M e>n (t + At,t) ) 
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For any pair of time (t, s ) € Fq, assume m = Kj G IN, we have 


m 

n Mn ( s +*At,«+(*—i) At) 

Z=1 

V (5 n (At)r 

\ o M e>n (t,s ) 

where S^Jt - h) = (5' n (At)) m - i S' Xjn (A*)S’j /!ri (Af)5 z , n (s + zAt - /i), h G [s + (z - 
1) At, s + zAt]. In order to show that n"=i M n (s + zAt, s + (z — 1) At )P n is strongly 
convergent to P n M(t, s), we only need to show 

ll'S'ce,n( i - h)M e>n (h, s)P n f - P n S{t - h)M e (h,s)f || ->■ 0, m,n ->• oo. 



In fact, 


||Sje, n (t - h)M e , n (h, s)P n f - P n S(t - h)M e (h, s)f\\ 

< ||sje,„(« - S)P„/ - Si e n (t - h)P n M e (h, s)/|| 

+ IISje,n(« - MPnM e (/z, s)/ - P n S(t - fc)M c (fc, «)/||, 

where, clearly, ||S'* en (t - h)M ejn (h, s)P n f - S l cen (t - h)P n M e (h, s)/|| ->• 0, m,n ->• 
oo. Moreover, we have 


||5^ in (t - h)P n M e [h, s)f - P n S(t - h)M e (h, s)/|| 

^||((5 n (At)) m -' — S(t — s — iAt))S x , n (^)S y , n (^)S z , n (s + iAt - h)P n M e (h, s)/|| 

+ || S(t - s - iAt)(S x , n (^)Sy jn (^)S Ztn (s + iAt - h)P n 
- P n S(s + iAt - h))M e (h , s)/||, 


where, according to (EB, || ((5„(At)) m 1 - P n S(t - s - iAt))S Xtn (^)S y , n (^£)S Z!n (s + 
iAt — h)P n M e {h , s)f || —> 0 and 

\\{Sx,n(^)Sy,n(^)S z ,n(s + zAt - h)P n - P n S(s + iAt - h))M e (h,s)f II 

^ \\S x ,n(^-) S y,n(^-)S z ,n(s + iAt - A) II 

•11(1 - S z , n (h -s-(i- l)At))P n M e (h,s)f\\ 

+ \\(S x , n (^)Sy, n (^)S z ,n(At) - S n (At))P n M e (h, s)f\\ 

+ \\(S n (At)P n -P n S(At))M e (h,s)f\\ 

+ \\(S(h — s — (i — l)At) — I)S(s + iAt — h)M e (h, s)/|| —> 0, At —» 0, n —> oo. 


Further, we discretize the model in time by the Lax-Wendroff scheme for advection 
equations in horizontal directions and Crank-Nicolson scheme for the diffusion equation in 
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the vertical direction such that S x / y / Zjn are approximated by 


°x/y,n\ 0 ) 1 ' 0 11 x/y,n < o A x/y,n^ 


2 
At 
~2 

At 


Sz,n{At) = (I- ^D z , n )-\I + ^-D z 


B e z , n {t,8)f = (I-=fD Ztn 


)-\^-M e , n (t,s)f), 


where A x j yn and I) z , n is the approximate generators to n-dimensional state space based 
on finite difference methods. 

It is well known ifTTft that the Lax-Wendroff scheme is consistent and conditional sta¬ 
ble for A x and A y and the Crank-Nicolson scheme is consistent and stable for D z , (I — 
is the consistent and condition stable implicit Euler scheme, by Lax equiva¬ 
lence theorem, that is 


lim Q \\(Sx/y/z,n( At )) At f ~ Sx/y/z,n(t)f || =0, / G L 2 (^), 

lim ||((7 - ^ D 2 , n )- 1 )^/ - S x , n (t)f\\ =0, / G L 2 (Q). 

\t-> 0 Z 


Similarly defining S n .— Sxin$y,n$z,n$y,n$x,ni 


lim ||(S n (At))^P n /-P n S(t)/||=0, / G L 2 (D). (22) 

n—>-oo,Ai—>-0 


~ t 

Since A x , A y and D z are self-adjoint, S*(At)) Ai is also strongly convergent to S*(t). 
Thus, (fl9l ) is approximated by 

/ SCn(t + At) \ 
v Se n (t + At) ) 

( S n (At) S x , n (%)S y , n (%)Bl n (t + At,t) \ f Scn(t) \ 
v o M e>n (t + At,t) ) V 5e n (t) ) 

^ V(t)V(t)(^ " ^D^-^iddnit + At) + 5dn(t ))] ^ 


Defining the above block evolution operator as M n (t, s), (t, s) G Tq, we have 


m 

M n (s + iAt, s + (i - l)At) 

i= 1 

( (Sn(At)r 5Xi(5 B (At))"*-‘5* >n (^)5 yiB (^)fl5 >B ( a + iAt ja ) 
V 0 M e n (t, s) 


We define by 


B e zn (s + iAt , s + (i — 1) A t, s)f 
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L 


•s+iAt 


S z ,n(s + iAt - h)M e) n(h, s)fdh, f € L 2 (yi). 


J 1) At 

By the trapezoidal rule and convergence of the implicit Euler scheme, we have 




For the observation system, we assume there is only a single observation during the 
entire time interval and define the observation mapping H r : L 2 (Q) —> R by 



where f! r and V r arc similarly defined as (l20l) . Then, the observation system extended with 
the emission rate is given by 



where 5y(t) € E and u(t) is the white noise with distribution iV(0,1). 

According to the spatial discretization of the model, in the vertical direction, [0,1] is 
discretized into three layers {0,0.5,1}. Since the diffusion coefficient K(z) is small, we 
assume possible locations of the single observation are around the grid points in the first 
layer z = 0. 

We have already shown that the assumptions (Al) — (A3) in Section [6] and the com¬ 
pactness of the possible area of observation locations are satisfied. 

In addition, according to the spatial discretization, we assume that the initial covariance 
is given by P n (to\t-i) = e~ 8 I n , where I n is the n x n identity matrix. It implies that 
P n (to\t-i) does not converge to a nuclear operator. It is shown in Figure 1 that the optimal 
location and minimal cost based on Kalman filter do not converge in this situation. 

Then, according to (1221) and dominated convergence theorem, we obtain that 



i =1 


is strongly convergent to 



2=1 
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Minimal cost at t=3 by KF 


Optimal location at t=3 by KF 
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Figure 1: Minimal cost and optimal location based on Kalman filter without the nuclearity of P(fo|f-i 
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Figure 2: Minimal cost and optimal location of the estimation of the state at final time by Kalman filter. 
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Minimal cost at t=0 by KS 
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Figure 3: Minimal cost and optimal location of the estimation of the initial state by Kalman smoother. 
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Further, -=i M n (s + iAt, s + (i — 1) At) is strongly convergent to n”=i M n (s + iAt, s + 
(i - 1)A t). 

Next we define the initial covariance as 

OO 

P(to\t-i)f = J2e- i2 (f, ei ) ei , f € L\n), 

i =1 

where {e*} is an orthogonal basis of L 2 (fl). The n-dimensional approximation of P(to|t_i) 
is given by 

n 

P n (f 0 |f_i)P n / = Y, e ~ i2 ( Pn/, G>G, / € T 2 (Q). 

i —1 

With this choice, P(to|t_i) is nuclear and the assumption {AA) in Section[6]is satisfied. By 
Theorem 16.41 the optimal location and minimal cost based on Kalman filter and smoother 
are convergent, which are shown in Figure 2 and Figure 3, respectively. 


References 

[1] A. BATKAI, P. CSOMOS. AND G. NICKEL, Operator splitting and spatial approxi¬ 
mations for evolution equations , J. Evolution equations, vol. 9, pp. 613-636, 2009. 

[2] A. BATKAI, P. CSOMOS. B. FARKAS, AND G. Nickel, Operator splitting with 
spatial temporal discretization, Oper. Theory Adv. Appl. 221. pp. 161-171, 2012. 

[3] R. BUIZZA, Localization of optimal perturbations using a projection operator, Quart. 
J. Roy. Meteor. Soc., 120, pp. 1647-1681, 1994. 

[4] R. BUIZZA AND A. MONTANI, Targeting observations using singular vectors, J. At¬ 
mos. Sci, vol. 56, pp. 2965-2985, 1999. 

[5] D. E. CATLIN, Estimation, control, and the discrete Kalman filter, Springer-Verlag, 
1989. 

[6] R. F. CURTAIN and A. J. PRITCHARD, Infinite Dimensional Linear Systems Theory, 
Lecture Notes in Control and Information Sciences, Vol. 8, 1978. 

[7] R. F. CURTAIN and H. Zwart, An Introduction to Infinite-Dimensional Linear Sys¬ 
tems Theory, Berlin, Germany: Springer Verlag, 1995. 

[8] R. CURTAIN AND A. SASANE, Compactness and nuclearity of the hankel opera¬ 
tor and internal stability of infinite-dimensional state linear systems, Int. J. Control, 
vol. 74, no. 12, pp. 1260-1270, 2001. 

[9] R. Curtain, K. Mikkola, and A. Sasane, The Hilbert-Schmidt property of feed¬ 
back operators, J. Math. Anal. Appl., vol. 329, pp. 1145-1160, 2007. 

[10] N. Darivandi, K. Morris and A. Khajepour, An algorithm for LQ optimal 
actuator location, Smart Mater. Struct. 22, Jan. 28, 2013. 


30 


[11] R. DAUTRAY, J. L. LIONS, Mathematical analysis and numerical methods for science 
and technology: Volume 6 Evolution Problem 2, Springer, 1999. 

[12] H. Elbern, A. Strunk, H. Schmidt and O. Talagrand, Emission rate and 
chemical state estimation by 4-dimension variational inversion , Atmos. Chem. Phys., 
vol. 7, pp. 3749-3769, 2007. 

[13] P. L.Falb, Infitite-dimensional filtering: The Kalman-Bucy filter in Hilbert space. 
Info, and Control, vol. 11, pp. 102-137, 1967. 

[14] J. C. GEROMEL, Convex analysis and global optimization of joint actuator location 
and control problems, IEEE Trans. Autom. Control, vol. 34, no. 7, pp. 711-720, 1989. 

[15] J. GIBSON, The Riccati integral equations for optimal control problems on Hilbert 
spaces, SIAM J. Control Optim., vol. 17, no. 4, pp. 637-665, 1979. 

[16] J. GIBSON, Linear-quadratic optimal control of hereditary differential systems: Infi¬ 
nite dimensional Riccati equations and numerical approximations, SIAM J. Control 
Optim., vol. 21, no. 1, pp. 95-139, 1983. 

[17] C. D. JOHNSON, Optimization of a certain quality of complete controllability and 
observability for linear dynamic systems, Trans. ASME, vol. 91, series D, pp. 228- 
238, 1969. 

[18] R. E. KALMAN, A new approach to linear filtering and prediction problems. Jour, of 
Basic Engineering, pp. 35-45, Mai - ., 1960. 

[19] R. E. KALMAN and R. S. BUCY, New results in linear filtering and prediction the¬ 
ory, Jour, of Basic Engineering, pp. 95-108, Mai - ., 1961. 

[20] E. N. LORENZ, A study of the predictability of a 28 variable atmospheric model, 
Tellus, vol. 17, pp. 321-333, 1965. 

[21] R. K. MEHRA, Optimization of measurement schedules and sensor designs for linear 
dynamic systems, IEEE Trans. Autom. Control, vol. AC-21, no. 1, pp. 55-64, 1976. 

[22] K. MORRIS, Linear-quadratic optimal actuator location, IEEE Trans. Autom. Con¬ 
trol, vol. 56, pp. 113-124, 2011. 

[23] C. Nam, Y. Kim AND T. A. Weishaar, Optimal sizing and placement of piezo- 
actuators for active flutter suppression, Smart Mater. Struct, vol. 5, pp. 216-224, 1996. 

[24] S. L. PADULA AND R. K. KINCAID, Optimization strategies for sensor and actuator 
placement. Technical Report NASA/TM-1999-209126, Apr., 1999. 

[25] K. R. PARTHASARATHY, Probability measures on metric spaces, Academic press, 
1967. 

[26] J. Weidmann, Linear Operators in Hilbert Spaces, New York: Springer-Verlag, 
1980. 


31 



[27] N. N. Yanenko, The method of fractional steps: solution of problems of mathematical 
physics in several variables , Springer, 1971. 


32 



